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We formulate a model for a cooperative ballistic deposition (CBD) process whereby the incoming 
particles are correlated with the ones already adsorbed via attractive force. The strength of the 
correlation is controlled by a tunable parameter a that interpolates the classical car parking problem 
at a — 0, the ballistic deposition at a — 1 and the CBD model at a > 1. The effects of the correlation 
in the CBD model are as follows. The jamming coverage q(a) increases with the strength of attraction 
o due to an ever increasing tendency of cluster formation. The system almost reaches the closest 
packing structure as a — > oo but never forms a percolating cluster which is typical to ID system. In 
the large a regime, the mean cluster size k increases as a 1 / 2 . Furthermore, the asymptotic approach 
towards the closest packing is purely algebraic both with a as 5(00) — q(a) ~ a~ 1//2 and with k as 
5(00) — q(k) ~ fc" 1 where q(oo) ~ 1. 

PACS numbers: 05.20.Dd,02.50.-r,68.43.-h 



I. INTRODUCTION 

The kinetics of a monolayer growth by the deposition 
of macromolecules and colloidal particles onto solid sub- 
strates has been the subject of extensive research for the 
past years (see jj], ||, U for extensive review) . The reason 
is well justified because its importance and significance 
cover many seemingly unrelated topic in physics, chem- 
istry, biology and other branches of science and technol- 
ogy. From a theoretical point of view, the random se- 
quential adsorption (RSA) of a monodisperse particle is 
one of the simplest model that can describe deposition 
phenomena Q. In this process, particles are deposited 
randomly, one at each time step, with the strict restric- 
tion that overlapping is forbidden. This can be described 
by the following algorithm, (i) At each time step, a ran- 
dom position is chosen from the whole substrate and is 
assigned to the center of the particle picked for deposi- 
tion, (ii) If the incoming particle collides with a previ- 
ously adsorbed one, the trial attempt is rejected; other- 
wise it is adsorbed irreversibly, (iii) In either case, the 
time is increased by one unit and the steps (i) - (ii) are 
repeated until the system reaches a state when particles 
can no longer be adsorbed. 

One of the virtue of the RSA model is that like many 
statistical physics problems, it is exactly solvable in one 
dimension in both its continuum and lattice version for 
some specific cases. The one dimensional continuum ver- 
sion of the model is popularly known as the random car 
parking (RCP) problem and has attracted much atten- 
tion. Despite the inherent simplicity in the RSA model, 
it still captures essential generic features of the process 
and has proved to describe successfully the behavior of 
many experimental systems, namely, the adsorption of 
proteins, latex and colloidal particles Jj| [| [?J . Neverthe- 
less, there have been continuous research efforts to in- 
clude various important physical features to make it more 
realistic and thus covering a wider range of real life situ- 
ations H |9[ [l(J. Along this road, a good deal of progress 



has already been achieved and yet we are far away from 
a complete theory. In recent years it has received extra 
momentum and the number of papers published in the re- 
cent years is a clear testimony to it || O, [l2| [ll| |l4|, [l5| . 

The strongest criticism of the RSA model is its out- 
right rejection of particles that fall upon an already ad- 
sorbed one. The most successful model overcoming this 
criticism is known as the ballistic deposition (BD) pro- 
cess Jll], 0. The BD model is best explained in terms 
of the deposition of disks of diameter m on a line. In 
this case, whenever an incoming disk overlaps with an 
already adsorbed one, it is allowed to roll over the lat- 
ter disk following the path of the steepest descent. In 
doing so, the disk can either touch the adsorbing plane 
(global minimum) or it may find itself trapped in the lo- 
cal minimum formed by two or more connected disks. In 
the former case, the disk is irreversibly attached with the 
one it rolled over leaving no gaps in between, while in the 
latter case the trial attempt is rejected. Both in the sim- 
ple RSA and the BD model, only a short range hard-core 
repulsion via the excluded volume effects is taken into 
account. All forms of long-range interactions between 
the particles in the adsorbed and adsorbing phases are 
ignored. There are some fragmented attempts though, 
to include some specific forms of interactions such as the 
electrostatic, dipolar and the hydrodynamic interactions 




In this article, we consider a model that includes the 
attractive force between the particles in the adsorbing 
and the adsorbed phase mimicking the long range in- 
teraction. To study the underlying mechanisms in such 
complex phenomena, like the deposition processes, it is 
of great advantage to have a flexible model which al- 
lows explicit analytical solutions with different possible 
mechanisms. In order to increase the flexibility of the 
model, we further generalise it by introducing a parame- 
ter a that can tune the strength of the attractive forces. 
This would certainly facilitate the study of the general 
effect of long-range interactions in the whole process. We 
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can recover the simple RSA results by setting a = 0, the 
BD results by setting a = 1 and while a > 1, we show 
that the model describes the cooperative ballistic depo- 
sition (CBD). Obviously, the strength of correlation is 
determined by the strength of the attractive force. This 
results in an increased overlapping with the already ad- 
sorbed particles. However, once a particle collides with 
an already adsorbed one, it follows the rule of the sim- 
ple BD. Thus, as the strength of interaction increases, 
we expect an increasing rate of successful adsorptions. 
The analytical part of the model we study in this work 
is similar to the one previously studied by Viot et al pl| . 
However, we go far beyond that not only by offering a 
completely new model but also by giving an in-depth 
analysis and an exact physical explanation as well as 
extracting more insightful results from several aspects. 
First, we give a general explanation of the model and a 
means to convert it into a simple and well known BD 
model based on a philosophical approach. Our ideas are 
backed up and well supported by direct numerical simu- 
lations which enables to understand the physical nature 
of the system described by each term of the rate equa- 
tions we used for analytical solutions. Second, we present 
more extensive results showing the asymptotic approach 
of the coverage towards the jamming limit in terms of 
the various parameters involved in the processes. 



II. CBD MODEL 

We consider a system which consists of a reservoir of 
particles with diameter to lying in the immediate vicinity 
of a substrate which is an infinitely extended line. The 
adsorbing particles may be in the gas or in the fluid phase 
and arrive in the adsorbing plane through Brownian mo- 
tion. As soon as a particle comes into contact with a 
gap large enough to accommodate it, it is then adsorbed 
immediately and irreversibly. In addition, the incoming 
particle that touches an already adsorbed one, is allowed 
to follow the BD rules to form a monolayer. The sim- 
plicity of the ID problem lies in the following situations. 
Every successful deposition of a particle on a given gap 
divides it into smaller gaps having the same geometry as 
the parent gap. It is this shielding property, found only 
in ID, which we shall use to gain further insight into the 
problem and tackle it analytically. For simplicity's sake, 
we assume that the daughter gaps are uncorrelated, irre- 
spective of the island size separating the gaps from their 
neighbours, so that we can treat each gap as an inde- 
pendent entity. We further assume that each roll-over 
motion is completed prior to the next trial attempt for 
deposition. 

At this point, it is useful to discuss the classical RSA 
and the simple BD model before introducing the coop- 
erative BD model. In the classical RSA, an incoming 
particle is adsorbed successfully in a gap of size x, if the 
center of the incoming particle arrives in any place but 
to/2 away from either edges of the gap. This means that 
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FIG. 1: Schematic illustration of the model in different situ- 
ations. The shaded regions between AB and CD in (a) have 
local minimums and play no role except for kinetic reasons. 
We therefore eliminate all such shaded regions so that we have 
a system where all the gaps are separated from their neigh- 
bours by only one disk as shown in (b). The dotted lines in 
(b) represent the probability distribution in different regimes 
of a given gap. The steepest descent path in (c) is artificially 
increased by using the idea of a virtual disk at the expense of 
lowering the height of PD so that every point of x + (2a — l)m 
is now equally likely to be chosen by the trial attempt. 



only (x — to) of a given gap x is accessible for adsorption, 
which we have illustrated in Fig. 1(a). In the BD model, 
on the other hand, those particles that fall on an already 
adsorbed one may reach the substrate via the roll over 
motion. The deposition via rolling is successful if the 
center of the incoming particle falls within a distance of 
to/2 on either sides of both edges. It is then adsorbed on 
the respective edge creating a new gap of size (x — to). 
That is, for a given gap x, the total position accessible 
to a new arrival is (x + to) which is shown in Fig. 1 (a) . 
Note that any particle dropping in the shaded regimes 
AB or CD are considered to be trapped due to the lo- 
cal minimum and will never reach the global minimum or 
the adsorbing substrate. Thus, if we are not interested in 
the kinetic aspect of the process, we can safely delete the 
shaded regimes as if they did not exist and assume that 
the neighbouring gaps are separated by only one disk as 
shown in Fig. 1(b). We can thus define each gap as an 
independent isolated interval bordered on either end by 
a semidisk so that if we had connected the two remote 
ends it would then form a ring with one particle at the 
joint. 

We are now in a position to introduce the long range 
attractive force among the incoming particles and the 
particles in the adsorbed phase. The question is how to 
incorporate it? First, we need to understand the effects 
of such an attractive force. The most significant one is 
that each adsorbed particle will tend to attract the in- 
coming particle towards it. This immediately breaks the 
random nature of the process as the incoming particles 
are more likely to land on an already adsorbed particle 
than on the gap. That is, as far as the selection of the 
trial position is concerned, we have two distinct proba- 
bility distribution (PD) regimes as indicated in Fig 1(b). 
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First, the force free regime where the attractive force is 
strictly zero. Second, the force field regime where the 
adsorbed particle exerts force and hence can influence 
the selection of the trial position. However, for analyti- 
cal simplicity, we assume a square well potential around 
the center of each adsorbed particle of width 2m. That 
is, the strength of the attractive force is equal anywhere 
within the force field regime. This is indicated by the 
dotted lines in Fig. 1 (b) representing the PD of the po- 
sition chosen by the trial attempt. The flat PD implies 
that all the points within that regime have the same a 
priori probability. 
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FIG. 2: Schematic illustration of the rules for the cooperative 
ballistic deposition model. The incoming disk is adsorbed 
directly as shown in (i) if it is dropped within (x — in). The 
deposition via the rolling mechanism is depicted in (ii). Here, 
any disk that falls on the steepest descent path of the virtual 
disk, OP, is assumed to be dropping on its corresponding 
equivalent point of the real particle QP from where it can 
successfully reach the global minimum. 

It is noteworthy to mention that the simple BD model 
refers to the case where both regimes have the same flat 
PD height and hence the whole substrate represents the 
zero force field or zero attractive force. Now, as soon 
as we switch on the attractive force, the height of the 
PD around each adsorbed particle will increase to a de- 
gree depending on the strength of the attractive force. 
However, we can lower the PD height by increasing the 
width keeping the total area unchanged. In this way, we 
can make the whole system having the same PD height 
and treat it like the simple BD model where an incoming 
particle can land anywhere in the substrate with equal 
probability including the disks that are already occupied 
as represented in Fig. 1(c). That is, we can envision it 
as follows: Prior to selecting a position for an incoming 
particle, we replace the size of the depositing particles 
used for defining the gaps by a virtual disk of diame- 
ter R > to, without altering the gap size. In this way, 
we increase the cross-section of collision for an incoming 
particle in order to count the additional collision proba- 
bility due to the attractive force. Once the position for 
the next trial deposition is chosen, we can immediately 
return to the system with disks having diameter m and 
proceed according to the simple BD rules as depicted in 
Fig. 2. In doing so, we artificially increase the proba- 
bility with which an incoming particle may collide with 
an already adsorbed disk thus enhancing the probability 
of adsorption via the rolling mechanism and mimicking 
the effect of the attractive force. One can thus expect an 



enhanced adsorption probability near the two extreme 
ends of each gap as the virtual disk size increases and, 
in the limit R — > oo, we can only expect the adsorption 
via the rolling mechanism except in the very early stage 
where the virtual diameter and the gap size may be of 
the same order in size. 



III. ANALYTICAL SOLUTION 

To address the problem described above analytically, 
we adopt the well studied rate equation approach of the 
gap size distribution function or concentration c(x,t). 
Let us define c(x, t)dx as the number of gaps at time 
t in the size range between x and x + dx. The kinetics 
of adsorption of the monodisperse particles can then be 
described by the following set of rate equations: 



dc(x, t) 
dt 



= -(x - m)c(x,t) + 2 / c(y,t)dy 

J x-\-m 

+2ara jc(ir + m, t) — c(x, £) j 



(i) 



for x > m and 



dc(x, t) 
dt 



J x+r 



c(y,t)dy + 2amc(x + m,t) (2) 



for x < to. The above rate equations are mean-field in 
nature as the fluctuations and correlations are ignored. 
The rate equation approach is based on the assumption 
that the creation and annihilation of gaps are indepen- 
dent of the size of the neighboring gaps. The first two 
terms of Eq. (1) and the first term of Eq. (2) are the 
same as that of the simple RS A process and thus describe 
the creation and destruction of a gap of size x due to the 
direct adsorption of size to on size y > x + m or on size 
x, respectively. The remaining terms in both equations 
also describe the creation and destruction of gaps but due 
to the rolling motion over the steepest descent path to 
travel the maximum linear path am, where a is a dimen- 
sionless constant number that we can tune. The factor 
'2' in the integral terms accounts for the fact that any of 
the two new gaps created upon a direct deposition on the 
gap size y > x + m can be of size x; whereas the same 
factor in the remaining terms takes into account that a 
gap of size x can be created or destroyed from either end 
by adsorption. In order to understand the role of a, it is 
convenient to rewrite Eq. (1) as 



dc(x, t) 
dt 



= -{x + (2a - l)m}c(x,t) + 2 
+2aTOc(a; + to, t). 



x+m 



c(y, t)dy 
(3) 



The term {x + (2a — l)m} in the above equation is the 
key to understand the role of a. Note that by setting 
a = we recover the classical RSA case where (x — m) 
of a given gap x is accessible for adsorption, which is 
consistent with our discussion in the previous section. 
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Similarly, a = 1 describes the simple BD model where 
the total positions accessible to a new arrival is (x + 
m) which is again consistent (see Fig. 1). In general, 
{x + (2a — l)m} means that a given gap of size x is 
bounded by a semidisk of diameter R = (2a — l)m while 
the adsorbing particles are of size m. That is, any particle 
that falls within a distance am from either end of {x + 
(2a — l)m} effectively will collide with the virtual disk. 
Every point of the steepest descent of the virtual path 
OP in Fig. 2 has its corresponding equivalent point on 
the real path QP. Therefore, an incoming particle falling 
on the virtual path is assumed as if it were falling on the 
exact equivalent position of the real path and vice versa. 

To solve Eq. (1) we seek a trial solution of the following 
form 

c(x,t) = A(t)e-( x - m ^\ (4) 

where A(t) is still an undetermined quantity fixed by the 
initial condition. Let us assume a monodisperse initial 
condition c(x, 0) = 5(x — L)/L so that we have 

lim / c(x, 0)dx = 0, lim / xc(x, t)dx — 1. (5) 

L — >ooj t >0j Q 

Substituting the trial solution into Eq. (1), we obtain 
the following differential equation for A(t) 

dlnA(t) 2e- mt n mt 

--^ = + 2ae mt . (6) 

dt t w 

Solving it, satisfying the initial conditions, we get 

A(t) = t 2 F(a,mt) (7) 

where the auxiliary function F(a, mt) is defined as 

F(a, mt) = e- 2 -lo mt i=^du+2a(i- ro t-e— *) _ (8) 

To obtain c(x, t) for x < m, we substitute the solution of 
Eq. (1) into Eq. (2) and then upon a direct integration 
we get 

c(x,t) = / u(2 + amu)F(a,mu)e~ xu du. (9) 
Jo 

The solutions c(x, t) can provide a complete analytical 
description of the process including its kinetic aspect. All 
we need now is to find useful ways of using these solutions 
for computing various physical quantities of interest such 
as the jamming coverage, the mean number density, the 
mean cluster size, etc. 

IV. NUMERICAL SIMULATIONS 

To test the physical description of our CBD model, 
we have simulated it on a computer. One obvious con- 
straint of the simulation is of course the finite size effect. 
However, for a sufficiently large substrate in comparison 
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FIG. 3: Schematic illustration of the recursive simulation 
scheme. 



to the depositing particles, the finite-size effect can be 
made sufficiently small. To simplify the simulation, we 
use the approximation mentioned in the previous section 
that the gaps are uncorrclated and can be treated in- 
dependently. This allows us to treat the problem in a 
recursive way if we focus on the coverage only and ignore 
the kinetics of the process. In brief, the description of 
the recursive scheme of the simulation is as follows. Let 
us assume that the initial gap of unit interval [0, 1] bor- 
dered on both ends by semidisks of radius (2a — l)m/2 
as shown in Fig. 3. We then generate a random number 
n from the interval [—(2a — I)m/2, 1 + (2a — l)m/2] to 
assign to the center of the incoming particle. The incom- 
ing particle is then adsorbed directly creating two new 
smaller intervals if n lies in the interval [m/2, 1 — m/2]. 
Otherwise, the disk is adsorbed on one of the edges of 
[0, 1] creating only one new gap of interval [m, 1] or of 
[0, 1 — m] depending on if is n € [—(2a — l)m/2, m/2] or 
n € [1 — m/2, 1 + (2a — l)m/2] respectively. In the case 
when the disk is adsorbed in one of the two edges, we 
increase the counter which provides the information on 
the mean cluster size of the system. We then continue 
the process assuming each new gap is again bordered by 
virtual semidisk and treat them in the same fashion as 
for the first step until we have no more gaps of size > m. 
At this point we add all the gaps of size < m and finally, 
using this, we can immediately calculate the jamming 
coverage and the mean cluster size. We have performed 
the simulation with substrate size ~ 10 6 m — 10 8 m within 
the interval [0, 1] and found an excellent match with the 
corresponding analytical results up to a several digits. 
We also noticed that increasing the substrate size by de- 
creasing the m value only contributes to a higher order 
precision as expected. 



V. RESULTS AND ANALYSIS 

The fraction of the line covered by the adsorbing par- 
ticles or the coverage 8(a, t) at different instants of time 
can be defined as 

/>oo 

9(a,t) = 1 - / xc(x, t)dx, (10) 
Jo 
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while the number density is defined by the following re- 
lation 



N(a,t) = / c(x,t)dx. 
Jo 



(11) 



However, we find it more convenient to handle their rate 
equation rather than their definition itself. The kinetic 
equation for the coverage is 



d9(a,t) 

dt 



= to / {x + m(2a - l)}c(x, t)dx — $(a,i). 



(12) 



Here, the quantity &(a,t) is the fraction of the substrate 
accessible to a new particle at a given time t. The kinetic 
equation for the number density on the other hand is 



dN(a,t) 
di 



(x — m)c(x, t)dx. 



(13) 



The above two equations can be combined together to 
obtain 



d8(a,t) 
dt 



dN(t) 
dt 



2am 



c(x, t)dx. 



(14) 



Note that in the simple RSA, one gap corresponds to one 
particle and thus we have 0(0, t) — mN(Q, t) reflecting 
the fact that the average particle size is the same as the 
size of the adsorbing particles. However, in the present 
case, for a > 0, the second term of the above equation 
describes the cluster formation. That is, the first term 
on the right hand side of Eq. (13) takes into account the 
direct deposition while the effect of the rolling mechanism 
is described by the second term. Using the solution for 
the appropriate boundary into Eq. ( |l2| ) yields 



0(a,t) 



F(a,u)(l + 2au)du. 



(15) 



This can provide all the information about kinetics aspect 
of the process, namely how the coverage evolves in time. 

One of the characteristics of the deposition process is 
that the system reaches a state of dead-lock in a finite 
time when particles can no longer be adsorbed. This 
is typically known as the jamming limit and the exact 
critical time to reach such a state should depend on a. 
However, till to-date there do not exist any theoretical 
means to pin down the exact critical time for reaching 
the jamming limit. Nevertheless, we can safely calculate 
the coverage in the jamming limit as 



q(a) = lim 0(a, t). 



(16) 



The jamming coverage has been of special interest in the 
study of the deposition phenomena as it can uniquely 
characterize the structure of the resulting monolayer. 
From the exact expression for the coverage, Eq. (15), 
it is of particular interest to know the approach of the 
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FIG. 4: Jamming coverage as a function of a that measures 
the strength of attraction: q(a) vs a. 




FIG. 5: The linear fits of ln(g(oo) — q(a)) vs ln(a) having slope 
— 1/2 in the large a regime reveals that the convergence of the 
jamming coverage towards the closest packing is power-law in 
nature. 



coverage 6(a,t) to the corresponding jamming limit q(a). 
We find that beyond the transient behavior, the system 
reaches its corresponding asymptotic coverage namely 
the jamming limit exponentially, with a decay constant 
2a, multiplied by an algebraic prefactor t^ 1 



q(a)-9(a,t) 



r x e- 2at , 



(17) 



which was also reported in pjj |. Obviously, the classical 
RSA (a = 0), we recover the power-law behavior which is 
also known as Feder's law p2| . Here for a > 0, the expo- 
nential approach towards the jamming limit reflects the 
fact that the increasing number of particles that land on 
an already adsorbed particles are successfully accommo- 
dated via rolling. Another interesting point to check is 
how the jamming limit varies as we increase the strength 
of interaction a. In other words, we want to see how 
the jamming limit changes as we increase the degree of 
correlation between the particles in the adsorbed and ad- 
sorbing phases. Fig. 4 shows a sharp rise in the jamming 
coverage at low a and a slow rise towards the closest pack- 
ing in the large a regime. In an attempt to quantify the 
slow regime we plot ln(g(oo) — q(a)) against ln(a) in Fig. 
5 and find that the jamming coverage converges towards 
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the closest packing obeying a power law 
q(oo) - q(a) ~ a~ 1/2 , 



(18) 



where q(oo) w 1. It is important to note here that the 
system never reaches a complete closest packing (q = 
1) even for a — ► oo. This is due to the fact that the 
substrate size too is of the same order as that of the 
size of the virtual disk, hence there is always a non-zero 
probability for a direct deposition at least in the early 
stage. We attempted to check it in the computer choosing 
both the virtual diameter and the initial substrate size to 
be of the same order and large enough to minimize the 
finite size effect. This is exactly the case described by the 
analytical model as we let a — ► oo. In doing so we never 
find a cluster covering the whole substrate . We checked 
it over and over again by increasing the substrate size 
and the virtual disk size up to ~ 10 1 0m. Nevertheless, 
neither analytical solution nor the simulation could give 
us an exact estimate for q(oo). 
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FIG. 6: Time dependence of the number density as a function 
of the coverage: The curves from top to bottom represents 
N(a, t) vs 6(a, t) for a = 0, 1, 5 and 10. 

We now intend to obtain an exact expression for the 
number density by substituting the solution c(x,t) for 
x > m into Eq. (H3) which yields 



N(a,t) = 



in™* F( a i u)du 



(19) 



The above relation for the number density immediately 
implies that it depends on the size of the adsorbing par- 
ticles. Here it is noteworthy to mention that the mean 
number density for the classical RSA (a = 0) is simply 
the coverage divided by the size of the adsorbing parti- 
cles. In this case, the mean number density increases lin- 
early with the coverage as time proceeds having a slope 
m during the process (see Fig. 6). However, as soon 
as a particle that fell on a previously adsorbed one can 
roll over the latter, the linear relation between N(a, t) 
and 9{a, t) is immediately ceased and it is replaced by 
a non-linear relation (see Fig. 6). Therefore, the mere 
knowledge of one of the two is not sufficient to obtain the 
other. This is due to the fact that the mean cluster size is 
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FIG. 7: The linear fits of the plot of ln(fe) vs ln(a) in the 
large a regime with slope 1/2 shows that the mean cluster 
size increases as ~ a 1 ^ 2 . 



different from the size of the adsorbing particle as the sys- 
tem keeps producing connected clusters of different sizes 
depending on the value of a. Fig. 3 shows that for a > 
the number density grows linearly at a very low coverage. 
This is due to the fact that, at an initial stage the incom- 
ing particles hardly encounter any pre-occupied species 
and therefore there exists almost no cluster. However, 
as the substrate gets crowded, it is evident from Fig. 6 
that the mean number density increases in a nonlinear 
fashion and the strength of non-linearity increases with 
increasing a. Therefore, to obtain the coverage (the num- 
ber density) from the number density (coverage) we need 
to know the mean cluster size. The expression for the 
coverage 9{a, t) and the mean number density N(a, t) at 
different instants of time can give us an estimate of how 
the mean cluster size defined as 



s(t) 



6(a,t) 
N(a,t) 



(20) 



grows in time and with the strength of a. We find that 
for a given size of the adsorbing species, the mean cluster 
size in the jamming limit is 



lim 



6(a,t) C F{a,u){l + 2au)du 

— S — — 7 

N(a,t) f™F(a,u)du 



(21) 



Obviously, like the mean number density, the mean 
cluster size should depend on the size of the adsorbing 
particle. However, the ratio between the two k = s/m 
remains constant in the jamming limit and therefore it 
is useful to call it the universal mean cluster size. It is 
note worthy to mention here that for the simple RSA we 
get k = 1 and hence the mean cluster size is the same 
as that of the size of the adsorbing species. However, for 
a > we find that the universal mean cluster k > 1 and k 
increases monotonously with increasing a. Furthermore, 
Fig. 7 reveals that in the large a regime the universal 
mean cluster size k increases as 



,1/2 



(22) 
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The probability of adsorption of particles without over- 
lapping with any preadsorbed particle decreases with 
time and for all a; however the strength of decrease gets 
sharper and sharper as a increases. For a — > oo, only 
at a very initial stage some particles may be adsorbed 
by direct deposition without overlapping. Thus, the sys- 
tem never reaches a state of closest packing but of almost 
closest packing. This is well supported by our numerical 
simulation. The jamming coverage thus increases with 
increasing k (see Fig. 8). We find that like the q{a) vs a, 
the approach of the jamming limit towards almost clos- 
est packing against k also follows a power law form but 
with a different exponent. As shown in Fig. 9, the plot 
of q(k(a — ► oo)) — q(k) versus k in the logarithmic scale 
along both axis is well fitted by a straight-line with slope 
1 and hence 

q(k(a — > oo)) -q(k) ~ k~\ (23) 
where q(k(a — ► oo)) w 1. 

1 I 1 1 1 T7T71 
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FIG. 8: The jamming coverage as a function of the mean 
cluster size k. 
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FIG. 9: The linear fits of ln(g(oo) — q(a)) vs ln(fc) in the 
large k regime shows that the jamming coverage increases in 
a power-law fashion ~ fe . 



VI. DISCUSSION 

We first discuss various possible interpretations for dif- 
ferent regimes constituted by different values of a in our 
CBD model. It is interesting to note that when a = 1, the 
size of the virtual and the real disk coincides and hence 
every point along the substrate has the equal probability 
of occupation by the incoming particle. In this case, par- 
ticles that fall on an already adsorbed one may travel up 
to a linear distance of its own size m via the rolling mo- 
tion following the steepest descent path. The particle can 
either be trapped in the local minimum or it can reach the 
global minimum. In the former case, it is rejected while 
in the latter case it is adsorbed irreversibly touching the 
one it just rolled over, provided there is at least a gap to 
accommodate it. The situation for < a < 1 is also in- 
teresting as it describes a mixed process composed of the 
simple RSA and the BD process. In this case, whenever 
an incoming disk encounters another disk, it is then al- 
lowed to roll over the latter disk with probability a. The 
rolling mechanism over the previously adsorbed disk can 
be assumed to be due to the gravitational pull towards 
the adsorbing surface. The ballistic deposition model 
(0 < a < 1) can thus describe the situation where the 
transport of the adsorbing species is dominated by gravi- 
tational effects. However, the model with 1 < a < oo can 
describe the situation where the gravitational effects as 
well as the attractive interaction between the elements of 
the adsorbed phase and the incoming particles play the 
dominant role. 

It is noteworthy to mention that Viot et al also gen- 
eralised the ballistic deposition model and gave the fol- 
lowing definition plL The disks of unit diameter are 
dropped uniformly and sequentially one at each time 
step. The disks can either reach the adsorbing plane 
or fall on an already adsorbed disk or on a cluster of 
disks. In the former case, the trial attempt is retained 
with a probability q = (1 — p). In the latter case, on the 
other hand, the trial disk follows the path of the steepest 
descent over the disk it encountered. The disk is then ad- 
sorbed with probability p provided the particle can reach 
the global minimum by the roll over motion; otherwise it 
is trapped in an elevated position and it is rejected. In 
pl[ the parameter a is defined as a — p/q. Therefore, 
the direct adsorption and the deposition via the rolling 
mechanism will both cause an increase in the rejection 
rate of the trial attempts for deposition. As a result, 
the approach to the jamming coverage in time should be 
slower as a increases. However, the analytical solutions 
reveal an opposite behaviour, see Eq. (17). 

Note that it is the rejection criterion that determines 
how fast or how slow the coverage should evolve and fi- 
nally reach the jamming limit. We would like to point 
out here that the apparent ambiguity arises because one 
cannot define a = p/q due to the following reasons. We 
can rewrite Eqs. (1) and (2) for a = p/(l — p) upon mul- 
tiplying both equations by (1— p). If we then set p = 1 in 
these equations we find that the remaining terms alone 



8 



are incapable of describing any meaningful physical pro- 
cess. The point to emphasize here is that the roll over 
mechanism comes after the trial attempt for deposition 
is made. The trial attempt to deposit the particle is the 
primary event at each time step of the process which may 
then be followed by the roll over motion upon collision. 
Therefore, the two events, the direct adsorption and the 
deposition via rolling are not mutually exclusive. 

The model we have presented in this article is solved 
exactly by means of a kinetic equation approach and sup- 
ported by the numerical simulations based on the recur- 
sive algorithm. The exact match between the analytical 
solution and the numerical simulation helped not only 
to confirm the validity of the mean-field approximation 
but also to shed a deeper insight into the nature of the 
problem. The basic principle of the model is the same as 
that of the simple ballistic deposition process. We have 
extended the simple BD by adding a certain degree of 
correlation between the adsorbing particles and those al- 
ready adsorbed. To increase the flexibility of the model, 
we allowed a parameter a that can tune the strength of 
the correlation which is induced by the attractive force. 
Instead of using the attractive force directly, we have 
shown a way of transforming it into a virtual situation, 
which is then just the simple BD model. 

The most significant consequence of the presence of the 
attractive force is that it results in an increased packing 
fraction and the mean clsuter size due to the formation 
of the higher order connected clusters. Moreover, the 
jamming coverage increases with the increasing degree of 
strength of the attractive force. Similar results has also 
been recently reported by Pastor-Satorras and Rubi Q , 
who studied a model of correlated sequential adsorption 
by numerical simulation both in one and two dimensions. 
However, unlike the square well type PD studied here, 
they used a Gaussian and exponential type PD around 
the center of each adsorbed particle. Nevertheless, de- 
spite the apparent differences in the detailed nature of 
the forces or in the PD, the qualitative behaviour seems 
remarkably identical to what we have found in this arti- 
cle. Pastor-Satorras and Rubi too observed the similar 
trend of the increase of jamming coverage as well as the 
mean cluster size. In addition, they too reported the 
approach towards the closest packing in the limit where 
the correlation is maximal. This reveals that the qual- 
itative effect of the attractive force is insensitive to the 
detailed nature of the force. However, in addition to solv- 
ing the model analytically we were able to quantify the 



effect of the attractive force. To this end, we have shown 
that in the strong force regime the convergence towards 
the closest packing against the strength of the attractive 
force follows a power-law relation q(oo) — q(a) ~ t~ 1//2 . 
Perhpas, the emergence of such a power-law behaviour 
implies a universal nature of the phenomena including 
the exponent in the sense that it is independent of the 
detailed nature of the attractive force. However, at this 
point it is just a conjecture and we intend to investigate 
it in our future work. 

VII. CONCLUSION 

In this article, we have presented an extension of the 
simple BD model by incoporating an attractive force be- 
tweent the elements of the adsorbed phase and the in- 
coming particle. The most significant consequence of the 
presence of the attractive force is the increase in packing 
fraction or the jamming coverage of the resulting mono- 
layer as we increase the strength of the attractive force a. 
This is manifested through the increase in the mean clus- 
ter size k and in the strong force regimes it increases as 
~ a 1 / 2 . We have shown that the system exhibits a power- 
law approach of the jamming coverage towards the clos- 
est packing both in a and k but with different exponents 
~ a -1 / 2 and ~ k~ x respectively, except for the weak 
field. It is important to note that although the jamming 
coverage increases with the degree of correlation vis-a-vis 
the mean cluster size, we can never create one single con- 
nected cluster spanning the whole substrate and giving 
the coverage q = 1 which is indeed typical to ID problem. 
Nevertheless, it indicates the potential structural phase 
transition in heigher dimensions and it has been indeed 
observed in (23) . 

Finally, as our model produces connected clusters of 
different sizes, it may well explain some aspects of RSA 
of polydisperse mixture in some sense. However, the mix- 
ture of particles has the following restrictions. Namely, 
the mixture contains a strict lower and upper cut-off and 
all the particles of intermediate size are of the integral 
multiple of the smallest particle in the mixture. 
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